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DESIGN OF STRUCTURES FOR OPTIMUM GE OMETRY 

Garret N. VANDERPLAATS, Research Scientist 
NASA, Ames Research Center 
Moffett Field, California, USA 

1. Introduction 

Optimization of finite element structures has received con- 
siderable attention in recent years. The majority of this work has 
dealt with structures of specified configuration, the design varia- 
bles, the cross-sectional areas of bar elements and the thicknesses 
of membrane panels. Relatively little effort has been directed 
toward determination of the optimum configuration of the structure. 
However, that work which has been reported [1 - 10] is sufficient to 
demonstrate that major design improvements can be achieved by allow- 
ing for configuration changes during the automated design process. 

A general method is presented here for geometric optimization 
of finite element structures. It is assumed that a reasonable 
initial geometry is specified. The number of finite elements, the 
number of Joints and the element-joint relationships are specified 
and are not changed during the optimization process. T1 j structure 
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may be statically Indeterminate and may support multiple loading 
conditions. Design variables Include geometric and member sizing 
parameters. The design objective may be minimum weight or coat, and 
constraints Include strength and stiffness limitations. The optimi- 
zation procedure is a general mathematical programming approach. 

The general optimization problem is formulated and the itera- 
tive design algorithm is outlined. This formulation is specialized 
for truss structures and design examples are presented. Conclusions 
based on this study are also presented. 

2 . General Formulation 

The general design problem considered here is to minimize the 
weight or cost of the structure. Constraints on the design can 
include allowable stress, buckling, displacement, dynamic and aero- 
elastic response limitations. The geometric location of the Joints, 
along with the element sizing variables, such as thickness or cross- 
sectional area, will be changed in order to optimize the structure. 
The general probl— .. ' . stated mathematically as: 

Minimize F(X) (1) 

Subject to: 

Gj (X) < 0 J - 1, m (2) 

where 



F(X) is referred to as the objective function. The set of m 
constraints which are Imposed on the design are defined by (2). The 
vector of design variables, X, includes geometric design variables, 
Xq, and member sizing design variables, X M that will be changed dur- 
ing the optimization process. The geometric design variables may be 
the coordinates of the joints themselves or may be the coefficients 
of any functional relationship which describes the geometry. For 
example, the surface of an aircraft wing may be described by a poly- 
nomial representation in order to conveniently interface the required 
aerodynamic analysis with the structural analysis. In this case, the 
geometry design variables may be the coefficients of this polynomial. 
Similarly, the element member .lzing variables, 5^, may be the actual 
thicknesses, of membrane elemt.its, or they may be the coefficients 



of a polynomial representation of the akin thickness distribution of 

a wing. 

In developing a geometric optimization capability, it is 
desirable to devise a technique which (1) deals with as few d« sign 
variables as Is practicable at any point In the opt lmlzat to', process, 
(2) will reduce the ill-conditioning Introduced hv mixing member 
sizing and geometric design variables, and (3) take full advantage of 
state-of-the-art techniques In fixed-geometry design. 

The approach used here is to treat the geometric design param- 
eters as independent design variables. The member sizing parameters 
are treated as dependent variables which are determined as a sub- 
problem. This is essentially the same as the approach used In [6], 
but Is extended here, using the techniques of [11], to deal with 
general finite element structures sub]ect to generalized constraints. 
Beginning with an initial geometric design vector, X^, the design 
proceeds iteratively by the following relationship: 

xf l • *8 ♦ • c 5 2 <*> 

where q is the Iteration number and S* 1 Is the search direction which 

* 

is yet to be determined, is a scalar parameter determining the 

distance of travel in the design space. 

For each proposed geometric vector, X^, the structure Is opti- 
mized with respect to the member sizing variables, X M , by solving the 
following suboptimization problem: 

Minimize F(X M > (5) 

Subject to: 

< 0 j - 1, m (6) 

Equations (5) and (6) are simply the standard form of the fixed- 
geometry optimization problem, and can be solved using any one of a 
variety r*f available algorithms. 

It is now necessary to determine the search direction, Sq. 
Assume that for the initial geometry the structure has been optimized 
with respect to the member sizing variables, and that for this sub- 
optimum design there are £ active constraints: 

Cj/X) - 0 k - 1, l (7) 

In practice, constraint G^CX) will be defined as active if its 
value is near zero, since precise zero is seldom meaningful on a 
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digital computer. It Is new necessary to find a search direction, 

S c , so that by moving In this direction in the geometric design 
space, the objective function will he reduced. This direction may be 
found by solving the following subproblem: 

Minimize VF(X) • S (8) 


Subject to: 


where 



(9) 

( 10 ) 


(ID 


Equations (8 - 10) themselves define an optimization problem in which 
the design variables are the components of S. This has the same form 
as the direction finding problem In the method of feasible directions, 
and the details for Its solution are described in [12] and [13]. The 
portion of S Is now substituted Into (4) and a one-dimensional 
search on is performed to updace the geomet 

The geometric optimization problem can .< ' o summarized in the 
following algorithm: 

1. Choose an initial set of design variables, X. Set q ■ 0. 

2. For the current set of geometric variables, X$, optimize the 
structure with respect to the member sizing variables, X M , using 
(5) and (6). 

3. Determine the number of active constraints (7) and obtain the 
gradients of the objective and active constraints (11). 

4. Determine the search direction, , using (8 - 10). 

5. Perform a one-dimensional search in accordance with (4) to deter- 
mine the value of a* which will minimize the objective function 

U 

subject to the constraints on geometry. For each proposed 
geometry, X^, update the member sizes using step 2. 

6. Check to see if the design has converged to the optimum. If no 
Improvement has been achieved, terminate. If the design has been 


4 



Improved, increase q bv 1 and return to step 3. 

3 . Application to Truss Structur e h 

The design algorithm is demonstrated here with application to 
three-dimensional statically Indeterminate truss structures. The 
structure is assumed to be linearly elastic and is subject to mul- 
tiple loading conditions. The objective is to minimize the total 
weight of the structure. Constraints Include limits on member siz- 
ing, element stresses. Joint displacements and Euler buckling of the 
elements. The weight objective function is 

NF. 

“ • E »i *i l i • < 12 

i-1 


where is the material density, A. the :ross-sectlonal area and 
the length of member i. NE is the total number of elements in the 
structure. Constraints on the design are defined as follows: 


Side constraints. 

> 

N 

1 

> 

1* 

0 

(13) 

Stress, 

v 1 

■*4 

o 

1 D 
•H 

D 

0 

(14) 



0 

(15) 

Displacement , 

4 lkJ / *lkJ ” 1 - 

0 

(16) 

Euler buckling. 

°lj /; bl ' 1 i 

0 

(17) 


where 


M 

’ci* °ti 


lower bound on member area i, 
calculated stress in member 1 under load condition j, 
compressive and tensile stress allowable, respec- 
tively, in member i. 


o. . ■ stress at which Euler buckling occurs in member i, 
bl 

6^ - displacement at Joint i in direction k (k ■ 1, 2, 3) 
under load condition j, 

■ allowable displacement at Joint l in directions k under 
load condition J. 

The stress at which Euler buckling occurs is defined by the relation- 
ship: 


- K iVi 


'bi 


(18) 
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where fcj la Young's modulus for element number i and is a pre- 
determined conetant which depends on the cross section of the bar 

element . 

The fixed geometry design utilized the constrained (unction min- 
imization program described In [14], This program is based on 
Zoutendljk's method of feasible directions (12) with modifications to 
deal with Initially Infeasible designs (11). Recent developments In 
approximation concepts by Srhmlt et al. [15 - 17] were utilized to 
gain maximum efficiency In the suboptimization problem. 

Three basic techniques were Incorporated here for sub- 
optimization. The first Is referred to as constraint deletion, where 
at each stage in the fixed geometry optimization problem only those 
active or near-active constraints are Included in the computations. 

The second Is the use of the reciprocal design variables to transform 
the constraint functions Into approximately linear form. The objec- 
tive Is nonlinear In reciprocal space, but is still explicit and is 
easily evaluated. Because of this approximate llnearltv of the con- 
straints, It Is logical to Include as the third technique a Taylor 
series expansion on the constraint functions which have not been 
deleted. This yields a linearized explicit form of the constraints. 
However, this is complicated somewhat by dealing with Euler buckling 
constraints as defined by (17). The stress at which Euler buckling 
occurs is a function of the design variables yielding a nonlinear 
constraint, even in reciprocal space. In this case, the Taylor series 
expansion is performed on the stress in member 1 under load condition 
J to yield an explicit linear form of the stress. The Euler buckling 
stress can then be reevaluated for each new proposed design, so that 
(17) is now explicit, but nonlinear. In the new design variables. 

The fixed geometry design proceeds by first analyzing the struc- 
ture and determining which constraints are active or near active. 

The appropriate Taylor series expansion is then performed on these 
constraints and this explicit form of the problem is used to mini- 
mize the weight. If the structure Is statically determinate, the 
solution is now complete. In the case of indeterminate structures, 
the structure is reanalyzed and the process Is repeated until con- 
vergence Is obtained. Even in highly indeterminate structures, this 
iterative procedure seldom requires more than 10 analyses, and a 
near-optimum design is usually obtained with only 5. 
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For geometric optimization, the Joint coordinates ol the 
structure were chosen as the design variables. 

4 . Design Exampl es 

A comj«:r*i program using these techniques was written fur truss 
optimization. The finite element displacement method is used for 
structural analysis and all gradient Information Is computed dlrectlv 
In closed form. The programming language Is FORTRAN IV and the 
examples were run on a CDC 7600 computer. 

Case 1 — 25-bar apace tnaa uith atreae ind i t • *i atrainta. 

Consider the 25-bar truss shown In figure 1. This structure has been 
used elsewhere as an example for fixed geometry design [IS] and con- 
figuration optimization [6]. 

z 

I 


l 90S m) *■ 



Figure 1 - Twenty-Five-Bar Space T'uss 

The truss was required to support two load conditions and was designed 
subject to constraints on member stresses and on Euler buckling. The 
loading Is the same as that used In [6], 

•6 2 

A minimum allowable cross-sectional area of 6.45*10 
(0.01 in. 2 ) was specified. The allowable stresses for all members 
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were specified ss 27.58*10^ N/nr (40,000 psl) In both tension snd 
compression. Young's modulus wss taken as 6.89*10^ N/m' ( 1 0 ^ psl), 
and the material denr.-iy p - 2767 kg/m' (0.1 lb/in. ^). The members 
were considered to be tubular with the nominal dlameter-to-thlcknes* 
ratio of D/t ■ 100 to give an Euler buckling constant (18) of K ( - 
39.274 for all members. The truss was required to remain symmetric 
with respect to both the x-z and the y-z planes. The Independent 
coordinate variables were taken as x^, , z ^ , x^ and y^. The remain- 

ing coordinates were linked to these design variables to maintain 
symmetry. The coordinates of Joints 1 and 2 were held constant, and 
Joints 7 through 10 were required to lie In the x-z plane. Member 

areas were linked In the following groups: A,, A,-A c , A -A., A, - 

1 l j by 1U 

*11' A 12~ A 13* A 14~ A 17* A 18~ A 21 an< * A 22~ A 2^* T ^ ere were then a total 
of five Independent coordinate variables and eight Independent area 

variables . 

The resulting optimum geometry is given in Table 1. For the 
final design, member 1 was constrained by its allowable tensile 

TABLE I. - 25-BAR TRUSS GEOMETRY 


Joint 



Coordinates In m 

(in.) 


Initial 

Final - Case 1 

Final - Case 2 

4 

X 

- 0.952 

0.544 

0.831 


1 

(37.5) 

(21.4) 

(32.7) 


y 

- 0.952 

1.176 

1.519 



(37.5) 

(46.3) 

(59.8) 


z 

- 2.54 

2.53 

2.939 



(100.0) 

(99.7) 

(115.7) 

8 

X 

- 2.54 

0.366 

1.255 



(100.0) 

(14.4) 

(49.4) 


y 

- 2.54 

2.123 

3.576 



(100.0 

(83.6) 

(140.8) 

stress. 

All 

other members 

were constrained 

by their respective buck- 


ling stresses In at least one member of each group. The weight of 
the truss was reduced by 48X from an optimum of 104.3 kg (229.9 lb) 
for the initial geometry to 54.4 kg (119.9 lb) for the final geometry. 
A total of 8 Iterations (4) on geometry were required to reach the 
optimum. A graph of weight versus iteration Is given in figure 2. 

The optimization required 27 f lxed -geometry designs using 171 analy- 
ses and 16.5 CPU seconds of computer time. A near-optimum design of 








61.0 k|( (1)4.) lb) was achieved In 4 geometry iteration* requiring 1) 

f lxetl -geometry dealgna with 81 analyse*. 


(Of 

o » »„ NitiAi. tiiin ocowt *»» oe* m.m 

• • ^ V * * t-MANGf 0 GtviMt tMV OetiMoW 



L -* A . _ A A 

0 2 « « a (0 

Iff RATION NOMHIR 

Figure 2. - Weight Versua Geometry Iter it ion Number 

Caae 2 —' j-bar apaoe tmaa iH.th otrcaa, bu<3«linni v. 1 diaplaoe- 

ment oonstrainte. The atructure dealgned aa Case 1 was reoptimized 
here with the additional requirement that the Joint dlaplacementa not 
exceed 0.00889 m (0.35 in.) in any of the coordinate direction* 
under any load condition. All other deaign conditions were the same 
as Caae 1. The optimum design for the initial configuration weighed 
255.6 kg (563.5 lb) and was constrained by the displacement limits at 
Joints 1 and 2. The optimum configuration weighed 74.3 kg (163.8 lb) 
and was achieved in 10 geometry iterations. The resulting geometry 
is given in Table 1 and the iteration history in figure 2. For this 
configuration, members 1, 10 and 11 were minimum size. The remaining 
cross-sectional areas were constrained by the Euler buckling limita- 
tion in at least one member of each group, with the exception that no 
stress or buckling constraints were active on members 6-9. Also, 
displacement limits were active on Joints 1 and 2. The design 
required 35 fixed-geometry optimizations using 204 analyses. The 
computer time was 20.9 CPU seconds. A near-optimum geometry weigh- 
ing 79.7 kg (175.7 lb) was achieved in 3 geometry iterations requir- 
ing 14 fixed -geometry optimizations and 72 analyses. 

5. Conclusions 

A general design algorithm has been presented for the optimum 
geometry design of finite element structures where a reasonable 
initial geometry has been specified. 

The basic conclusions of this study are: 

(1) Structural geometry can efficiently be treated as a design 


9 


parameter and major weight reductlona can often be achieved as a 
result of geometric change*. 

(2) By considering two separate design spate*, geometry and mem- 
ber sizing, the 111 conditioning problem* usually associated with 
combining member sizing and coordinate variables are seldom encoun- 
tered. Also, this approach allows for the application oi very effi- 
cient fixed geometry optimization techniques currently available. 

(3) The technique maintains the generality of mathemat leal pro- 
gramming, thereby allowing for the simultaneous consideration of a 
wide variety of constraints on the design. 

(4) The method converges rapidly to a near-optimum geometry 
(see fig. 2). 

(5) The computational technique Is in no way limited to trusses 
and applies directly to general finite element structures. 
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